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Abstract 

We present an elastic model of B-form DNA as a stack of thin, rigid plates 
or base pairs that are not permitted to deform. The symmetry of DNA and 
the constraint of plate rigidity limit the number of bulk elastic constants 
contributing to a macroscopic elasticity theory of DNA to four. We derive 
an effective twist-stretch energy in terms of the macroscopic stretch e along 
and relative excess twist a about the DNA molecular axis. In addition to 
the bulk stretch and twist moduli found previously, we obtain a twist-stretch 
modulus with the following remarkable properties: 1) it vanishes when the 
radius of the helical curve following the geometric center of each plate is 
zero, 2) it vanishes with the elastic constant K23 that couples compression 
normal to the plates to a shear strain, if the plates are perpendicular to the 
molecular axis, and 3) it is nonzero if the plates are tilted relative to the 
molecular axis. This implies that a laminated helical structure carved out of 
an isotropic elastic medium will not twist in response to a stretching force, 
but an isotropic material will twist if it is bent into the shape of a helix. 

PACS: 87.15.-v, 87.10.+e, 87.15.By. 



1 ohern@lubensky. physics .upenn.edu 
2 kamien@dept . physics . upenn . edu 
3 1 om@lubensky. physics . upenn . edu 
4 nelson@dept . physics . upenn. edu 



1 



1 Introduction 



The elastic properties of DNA have become a focus of recent research [[T], 
H, El> H H H- I n particular, stretching experiments on single molecules of 
DNA now provide a direct probe of the bending, stretching, and twisting 
elasticity of DNAfT], |[ pL0|| . Understanding the elasticity of single molecules 



of DNA may be relevant in vivo; e.g., the recA protein which is responsible for 
homologous recombination in bacterial mitosis has been observed to stretch 
and twist DNA when it is bound to the DNA molecule||. 

Marko and Siggia recently modeled DNA as a thin, uniform rod with a 
linear bending elasticity and calculated the extension of the rod as a function 
of the applied stretching force||. This worm-like chain model was appropri- 
ate for stretching forces / <C 7, where 7 ~ Ak B T/R 2 « 200pN is the stretch 
modulus, A « 50 nm is the bend persistence length, and R = 1 nm is the 
radius of the molecule. Forces in this regime pull on the thermal fluctuations 
of the molecular backbone but do not pull on the internal structure of DNA. 
The worm-like chain model accurately predicts the extension of the molecule 
to within 10% up to approximately 10 pN(J7|. Above 10 pN|7|, (SJ there are 
no longer thermal effects and DNA is stretched elastically. A strain vari- 
able e must be introduced to describe increases in the molecular length after 
thermally induced contour length fluctuations are fully stretched |2|, [11], [7|, §]. 

DNA is also characterized by the degree of twist of its phosphate back- 
bones about the central axis. Relative twist excess or deficit from a state 
of preferred twist is parameterized by a signed twist variable a. The chi- 
ral asymmetry of DNA allows for a twist-stretch coupling between a and 
e; this coupling is clearly visible in the recent experiments of Strick, et. 
q/.|T0|, |4|, [E|. Theoretical work has focused on calculating the extension of the 
molecule as function of both the applied stretching force and the applied ex- 
cess twist 0, |5], |T3|]. However, in this paper we investigate a more microscopic 



origin of the twist-stretch coupling [H 



We introduce a simple elastic model of B-form DNA in which the molecule 
is viewed as a stack of thin, rigid plates that represent base pairs and are 
rotated and displaced relative to each other. (See Section 3 below for a review 
of the geometry of B-form DNA.) The centers of mass of the plates define a 
helical path around a straight central line, which we call the molecular axis. 
The perpendicular distance from the molecular axis to the helix is called 
the helix axis offset. Distortion energies in this model are controlled by 



2 



an underlying continuum elastic energy with those elastic constants allowed 
by the anisotropic symmetry of DNA. Our model of DNA as a collection 
of thin, rigid plates has four contributing elastic constants coupling strains 
with spatial variation normal to the plates. Of particular importance is the 
elastic constant K 2 3 coupling compression perpendicular to the plates to a 
shear displacement parallel to the plates. K 2 z vanishes if the phosphate 
backbones point in the same directionF] and also vanishes in the isotropic 
limit in which each plate has reflection and C4 or higher symmetries. 

We calculate the twist-stretch coupling in terms of its elastic constants 
and the small helix axis offset. We find that the twist-stretch coupling has 
terms linear and quadratic in the helix axis offset. The coefficient of the 
linear term is proportional to the elastic constant i^23- In the isotropic limit 
or when K 2 3 vanishes, the twist-stretch coupling is quadratic in the helix axis 
offset as we calculated earlier [T2j in related but simpler models. 

An important simplification of our model is the constraint that the plates 
comprising the DNA stack are rigid and undeformable. This is equivalent to 
setting to infinity all elastic constants coupling the strains within a plate. We 



thus neglect propeller, buckle, opening, etc. deformations! 16] of individual 



base pairs and focus instead on inter-base pair deformations described by 
relative rotations (roll, tilt, and twist) and translations (slide, shift, and rise) 
of the base pairs |T5|, pj. We believe this is a reasonable approximation for 
stretching forces in the range 1 pN < / < 10 pN and relative twist excesses 
a < 0.05. 



2 Description of the Model 

DNA is composed of base pairs connected by two oppositely-directed sugar- 
phosphate backbones that wrap in two helices around the central molecular 
axis to produce the major and minor grooves. A cross-section of DNA normal 
to the phosphate backbones is shown in Fig.|l|. To construct our model for 
DNA, we first imagine unwinding the helix so that the phosphate backbones 
describe two straight, parallel paths. This unwound structure is an elastic rod 
characterized by an elastic constant tensor Kyki with components constrained 
by the symmetry of the rod. Each cross-section of the rod looks like Fig. 

5 Each phosphate backbone has a particular orientation: either 3'-5' or 5'-3'. 
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[I]. Let the z-axis be parallel to the phosphate backbones, and let the x- 
axis pass through the midpoint and perpendicular to the line connecting the 
backbones. The only symmetry operation of this rod is a rotation through 
7r about the x-axis that causes the axes to transform as x — > x, y — > —y, 
and z — > —z. This symmetry allows for a maximum of 13 independent 
elastic constants. Our assumption of the rigid shape of each cross-section in 
the x — y plane eliminates from consideration all strains that require spatial 
derivatives with respect to x or y. Thus we need only consider the strains 
u zz , u zx , and u zy and the four elastic constants Kz% = K zzzz , K\\ = K zxzx /4, 
K 2 2 = K zyzy /4:, and K23 = K zyzz /2 permitted by symmetry that couple 
these strains. The elastic constant K23 is allowed because the two phosphate 
backbones are oppositely-directed. If they were undirected or pointed in the 
same direction, the reflection y — > —y would be a symmetry and K23 would 
be zero. If, in addition, a rotation by tt/2, which transforms x — > y and 
y — > —x, were a symmetry, there would be only two elastic constants K 33 
and K = K\\ = K.22- We will refer to this case with two elastic constants as 
the isotropic limit. We view the cross-sectional slices as rigid plates that are 
coupled elastically via the elastic constants K33, K\%, K22, and K23. 

We now imagine that chiral forces inherent to DNA distort the straight 
rod to a helical structure in which the centers of mass of the rigid plates 
describe a helical path about a straight helical axis. This helical structure 
can be produced by displacing and rotating neighboring plates by a constant 
amount. We assume that chiral energies leading to a ground-state helical 
structure are linear in strain (i.e surface terms) so that the energies of dis- 
tortions from the ground state are determined by the elasticity of the original 
untwisted rod, i.e., by the elastic constants K33, Kn, K22, and K23. In the 
following sections, we will derive the twist-stretch coupling in terms of these 
constants and the helix axis offset. 

3 Geometry of B-Form DNA 

We now consider three different ways of stacking our thin, rigid plates to cre- 
ate a helical model of the straight state of B-DNA. We first imagine simply 
twisting the plates about the long axis passing through the center of mass 
of the original rod to create shape I. In this shape, the geometric centers 
of each slice lie on top of one another, the slices are rotated about the long 
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axis relative to one another, and the slices remain perpendicular to the long 
axis. Shape II is obtained by first twisting the plates about the long axis and 
then removing material from the side of each plate opposite to the phosphate 
backbones. The line connecting the new geometric centers of each slice fol- 
lows a helix. The slices remain perpendicular to the helix axis about which 
the geometric centers of each plate rotate. (In the discussion below we use 
the terms helix axis and DNA molecular axis interchangeably.) Shape III is 
obtained by bending the long axis of the sliced rod into a helix. The curve r 
connecting the geometric centers of each plate again follows a helix, but now 
the plates are perpendicular to the local tangent to f instead of the helix 
axis. 

We note that in the unstressed state of B-form DNA the plates are roughly 
perpendicular to the molecular axis, and their geometric centers nearly co- 
incide with the molecular axis[|T7]]. Thus to study small deformations of 
B-DNA, it suffices to restrict our attention to the case where shapes II and 
III are small perturbations of shape I. In shape II, the amount of material 
removed from each slice and hence the radius of the helix is small and in 
shape III the bend of the rod axis away from the molecular axis is also small. 
Shapes II and III can be described as twisted stacks of plates with small helix 
axis offsets. If the collection of plates has a helix axis offset, the plates are 
not stacked with one directly on top of the other; instead, the shift vector 
has a component perpendicular to the molecular axis. We will see that in 
our model the helix axis offset is the origin of the twist-stretch coupling. 

4 Elasticity Theory 

In what follows we will develop a long-wavelength elasticity theory for a 
collection of stacked, rigid plates in terms of small deviations of the shifts 
and rotations away from their unstressed values. We will then eliminate the 
shift and rotation variables in favor of the stretch along the molecular axis e 
and the relative excess twist about the molecular axis a. This will allow us to 
find the bulk twist and stretch moduli in terms of four elastic constants and 
geometric properties of the molecule. More importantly, it will also enable 
us to calculate the twist-stretch coupling and determine how it scales with 
the small helix axis offset. 

As shown in Fig. [I] we inscribe on each plate a right-handed, orthonormal 
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triad e a (n) where a = 1, 2, 3 and n is a unitless parameter which labels the 
plate. In the following, Roman indices run over Cartesian space coordinates 
(x, y, z) and Greek indices label the plate-fixed coordinate system (1,2,3). e x 
and e 2 lie in the plane of the plate, and e 3 is perpendicular to the plate. The 
coordinates on each plate are labeled by the pair (rji,rj 2 ) which corresponds 
to the point on the plate ^e^n) + r] 2 e 2 (n). We must also specify the origin 
on each plate. We choose the origin to be the geometric center (i.e., center 
of mass) of the plate, but this choice is arbitrary and does not affect the 
elasticity theory, as we will discuss in the Appendix. We can now describe 
the trajectory in space of the plates and the associated plate-fixed triads. 
We allow two sets of parameters: a vector of rotation rates u and shifts S. 
We decompose Q and S in the plate-fixed basis e a (n) and assume that the 
components u a and 5 a are constants independent of n. 

Q = oo a e a (n) 

7 . . where a = 1,2, 3. (1) 

o = o a e a (n), 

The vector of rotation rates uj describes the rate of change of the orientation 
of neighboring plates, and thus 

de 

— = to x e a = -e Q/ 3 7 ^e 7 , (2) 

— # 

where e Q/ g 7 is the antisymmetric tensor. The second set of parameters 5 
describes the relative displacement of two neighboring plates. We choose S 
to be the relative displacement of the origin f(n) of the (^1,^2) coordinate 
system. We have 

dr(n) 



dn 



5 = e a 5 a . (3) 



Note that according to our definitions u is unitless and 5 has units of length. 
Finally, the position in space of the point (771,772) on plate n is simply 
x(r]i,r]2,n) = r(n) + r/iei(n) +r/ 2 e 2 (n). For constant oj a and 6 a , these equa- 
tions will, in general, describe a helical structure. We may solve (0) and (|J) 
to find: 



r» = e,(0){n^^+ l^+^p^Mn) 



sin( \Q 








u 2 



(4) 
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Table 1: Geometric Properties of Several Helices 
Shape 5 uj Helix Axis Offset Pitch/27r 

i (0AS3) (0,0,^3) S3M 

II (0,S 2 ,S 3 ) (0,0, u 3 ) S 2 5 3 /u 3 

III (0,0, S 3 ) {0,u 2 ,u 3 ) u 2 S 3 / ^fuf+uji ^3/(^2 + ^3) 



with uj 2 = u) a u) a . Examination of (§) shows that the geometric center moves 
on average in the direction ep(Q)u}p/\uj\ and traces out a helix with helix axis 
offset r = \uj x <5|/|u;| and pitch p = 2%5 ■ uj/uj 2 . The shift vector, rotation 
vector, helix axis offset, and pitch are listed in Table [I] for each of the three 
shapes of B-form DNA we are considering. 

In a deformed state, the displacement and rotation rates change to 5' = 
5 + A(n) and uj' = uj + Q(n) and define new positions x'ijji, rj 2 , n) for points 
on plate n. Our goal is to calculate the energy of this deformed state relative 
to the equilibrium helical state to second order in the small parameters A 
and Q, which can in general depend on n. We consider plates of thickness 
Sn and calculate the energy in the deformed state of each plate to order Sn, 
ignoring terms of order (Sn) 2 or higher. We then sum over all plates to obtain 
the total energy of the rod. Since the energy of each plate is proportional to 
Sn, the sum over all plates can be converted into an integral over n. 

The energy of the nth plate will depend on the values of A and Q at 
n and, to the order we consider, not on their derivatives with respect to n. 
The energy of each plate will be the same function of A and Q for every n, 
so we need only calculate the energy of a single reference plate. Our elastic 
theory provides us with distortion energies of the reference plate as a function 
of the strains Uij = \(diUj + djUi) = x,y,z), where u = x' — x is the 
displacement variable and x,y, and z are the Cartesian coordinates of the 
reference plate. By a suitable choice of orientation we may take ei = x, 
e 2 = y, and e 3 = z for the plate of interest. The free energy for a single plate 
with thickness Sz is (in units of kBT): 



plate 



\ rz+5z r 

J J dxdydz K ijkl u i: jU k i. (5) 



kuT 2 
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Our goal is to express this energy as a function of A and Q and change 
integration variables from x, y, and z to 771, 772, and n. The Jacobian relat- 
ing the two coordinate systems is found by taking n a and n derivatives of 
%(ViiV2, n ) an d remembering that e a ■ = 8 a p. 

dxdydz = \5-z + prj^dnid^dn. (6) 

To calculate the strain tensor Uij we must take the derivatives of u(ni, r]2, n) 
with respect to n a and n and then relate these coordinates to Cartesian 
coordinates x, y, and z fixed on the plate. Since the derivatives of Ui with 
respect to 7] a only include terms proportional to 8n, the only derivatives in 
the strain tensor contributing to the total free energy in (J)|) are those in the 
n-direction. We therefore only need to calculate du/dn. Using the equations 
of motion, @ and (|]), and the expression for x(ni,T]2,n), we have to linear 
order in Q a and A a 



du(r] 1 ,r]2,n) 
dn 



A a + e a/37 ^?7 7 + (6 a + eapyUpty) e' a - e a . (7) 



We can neglect the last terms in (|7j) since they also are proportional to Sn. 
(We note that deformations of the base pairs can be included if we retain 
these Sn contributions.) 

We choose the slab at n to have its internal triad point along a Cartesian 
coordinate system (x, y, z) fixed on the plate with the z-direction aligned 
with e 3 and the x- and y-directions aligned with ei and e 2 respectively. To 
leading order in Q a and A Q , the only relevant components of the strain tensor 
come from the derivatives 

d z Ui = -q^-q^ = (A; + e ijk fljr] k )/6 3 , (8) 

where dn/dz = 1/63 for the three shapes we are considering. In (H) the sums 
over a, (3, and 7 have been replaced by sums over i, j, and k because we 
chose the body- fixed coordinates to be along (x, y, z). We have now reduced 
the number of effective elastic constants from 13 to 4 since only z derivatives 
contribute to the free energy. The free energy of a single plate in terms of 
the strains A and Q and the coordinates rji, 772, and n is: 

7^7= / dn f d 2 n ^ 3 + € ^ u ^ \K ij (A i +e Um tlir] m )(A j +e jkn Q k r] n )}, 

(9) 



8 



where Kij is related to K^i by 

( \K zizj i = j = x or i = j = y 
= I \K zizj i = y, j = z (10) 
l_ K z i z j i j z. 

We see that our elasticity theory of DNA as a stack of thin, rigid plates 
has four contributing elastic constants, K xx , K yy , K yz , and K zz . Below we 
will refer to these elastic constants as Ku, K 2 2, K 23 , and -^33, respectively, 
since the plate-fixed axes rotate in space as we move along the helical stack. 
We note that an isotropic theory has only two contributing constants, K = 
Kn/2 = K22/2 and ^33. The total free energy for the collection of plates is 
obtained by adding up the contributions to the free energy from each plate. 
The total free energy will have the same form as (|9|) with the Cartesian 
indices on the strains (x, y, z) replaced by the plate-fixed indices (1, 2, 3). 



5 DNA Effective Free Energy 

We now sum up the contributions from all slices and do the integrals over 
771 and rj2 in (|9|) to get an effective theory for the DNA molecule. Before 
we actually do this, we will make some further simplifications and approxi- 
mations. We first use the fact that under the DNA symmetry operation 772 
switches sign and hence all moments odd in r/ 2 are zero: (772) = (j]™ri2) = 
(m = 1,2,3, ...) where (•) = A^ 1 J d 2 r] (■) and A is the cross sectional area 
in the 771 — 772 plane. We also locate the origin of each plate at its geometric 
center so that (771) = 0. 

In the unstressed configuration of B-form DNA the base pairs make an 



angle of close to ir/2 with the molecular axisflTfl. We therefore take u\p Iq ^ 



63 where $3 ~ 3.4 A is the axial rise of the base pairs. In addition, the 



geometric centers of the plates are roughly located on the molecular axis JL7 
we therefore also assume ^1,2/^3 1. We will characterize the unstressed 
straight state of B-DNA by two known macroscopic parameters (the axial rise 
#3 and the twist rate of the plates about the molecular axis uj^ = 2tt5^/ Iq ~ 
0.63 rad) and two unknown microscopic parameters (62 and co> 2 )- I n what 
follows we will assume that 62/63, U2/UJ3, and the second moments (?7i)//q 
and (?7 2 )//o are small and work to quadratic order in these quantities. 
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The effective free energy is given below in units of ksT. We neglect the 
u 2 contribution arising from the Jacobian in (Q) since it yields terms that are 
third order in u> 2 /u 3 , {t]\)/Iq, and {t^/Iq. Including these terms does not 
alter our findings below for the twist-stretch coupling: 

^ = {A/28*) J dn {k x1 A\ + K 22 A\ + K 33 A 2 3 + K 33 (r,l)n\ + 

KM)Sl\ + ( K Kn(vl) + K 22 (r ] l}^nl-2K 23 (r ] 2 1 }n 2 n 3 + 2K 23 A 2 A 3 y (11) 

We note from ([□]) that the elastic constant K 33 acts as a stretch modulus 
and a bend modulus. 



6 Twist-Stretch Coupling 

Motivated by recent stretching experiments on torsionally constrained single 



molecules of DNA|T0f, we rewrite (|TT| ) in terms of the relative overtwist a 



about the molecular axis and the relative stretch e along the molecular axis. 
To accomplish this we must first find e and a in terms of the rotation and 
shift variables and Aj. After we incorporate e and a into (jTT|) , we minimize 
over the remaining unconstrained variables. This gives us the twist-stretch 
energy in terms of the elastic constants Kij and the geometric parameters 
of the DNA helix. We can then identify the stretch modulus B, the twist 
modulus C, and the twist-stretch modulus D. We note that this twist- 
stretch energy neglects the effects of thermal fluctuations; these effects have 



been studied recently!! p 



We now derive the twist-stretch coupling for the three shapes shown in 
Table [I] by considering the general case uj 2 , 8 2 7^ and then neglecting uj 2 8 2 
cross terms. The unstressed configuration is therefore uj = (0,uj 2 ,uj 3 ) and 
5 = (0, 82, 8 3 ). Each of the three equilibrium shapes mentioned previously can 
be obtained from this parameterization, for instance, shape II is generated 
by taking 8 2 7^ and uo 2 = 0. We define the stretch along the molecular axis 
as the relative deviation in the extension L along the molecular axis from its 
unstressed value L , 

e = A_ L (12) 

To find e in terms of Aj and fij we must determine the extension L from 
(f|). We see that each base pair step increases the extension by an amount 
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5 ■ uj/\u\ in the direction a;/|a;|. The total extension for N steps in the 
distorted configuration is therefore 

L = N5'~ (13) 
\uj'\ 

where the primes refer to the distorted state, i.e. 5' = (A 1; 5 2 + A 2 , S 3 + A 3 ). 
We can now expand ([12]) to linear order in the fli and Aj and to quadratic 
order in the small parameters <5 2 / £3 and UO2/0J3: 

^ _ A3 ,62 _ UJ2\ ^2 A2 Cj|^3 

#3 ^3 ^3 ^3 <$3 ^3 ^3 ' 

Note that when the helix axis offset of the helical stack is zero (62 = 002 = 0), 
the stretch reduces to e = A 3 /<5 3 . We also note that when aj 2 = A 3 and fl 2 
are the only distortions contributing to the stretch. 

We must express the relative excess link a = (Lk/Lko) — 1 in terms of 
the fli, where the linking number Lk is the number of times the plates rotate 
about the molecular axis over the length of the stack. In its unstressed state, 
DNA has linking number Lk = L /l . (We are working in the force regime 
where the molecular axis is nearly straight, and hence link and relative twist 
about the molecular axis are synonymous.) We see from @ that it takes 
27r/|u;| steps to complete one revolution about the molecular axis. Therefore 
in N steps the top plate has rotated by N\u\/2ii relative to the bottom plate. 
The relative excess twist about the molecular axis is therefore 

\ui\ 

where u' = (fii, uj 2 + fi 2 , u) 3 + fl 3 ). Expanding fli~5|) to linear order in the 0; 
and to quadratic order in uj 2 /uj 3 , we find 

a = ^ + (l-4)-- (16) 

UJ 3 UJ 3 V UJ% ' U 3 

Note that when u; 2 = 0, the base pairs are perpendicular to the molecular 
axis and the relative excess twist is simply VL 3 /uj 3 . 



Solving ([14]) for A3 and (|l6j ) for ^3, we substitute into ( |ITD and minimize 



over Ai 2 and f2i 2 . This gives the effective twist-stretch energy per length 
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along the molecular axis: 



^M = £(*. + ^ + 2D4 (17 ) 

where N is the total number of plates, B is the stretch modulus, C is the 
twist modulus, and D is the twist-stretch modulus. To zeroth order in the 
small parameters 62/ S 3 and 0^2/^3 the stretch and twist moduli are: 



22 



B = K 33 -K 2 2 JK 
C = (coi/5l)(K n ( V 2 2 ) + ( V l)(K 22 -K 2 23 /K 33 )). (18) 

B, C, and D have dimension L~ 3 and therefore the coefficients of the quadratic 
dimensionless strains in (|TT| ) scale as a persistence length divided by the 
square of the helical pitch as found previously||. The twist-stretch cou- 
plings Di for the three shapes given in Table |I] are shown below to lowest 
order in the helix axis offsets rf. 

D I = 

Du - 1^33 



^ 2 V 


)- 


-^22' 


'S 3 






^22' 


' s 3 



As 

An « ^(K a -^)f + (KM + K a -K a )(f)', (19) 
^33 v ^22 y S 3 V (r)() /V (5 3 / 

where the helix axis offsets for each shape are r\ = 0, rn = #2, and rm ~ 
uj 2 5 3 /uj 3 . In performing the calculation of the twist-stretch moduli we have 
made two simplifications: 1) we have dropped the quadratic terms in the helix 
axis offset that are proportional to K23 and 2) we have not included the u>2 
contributions arising from the Jacobian in (|9|). This last simplification does 
not affect our conclusions concerning the twist-stretch coupling (TS) shown 
in Table ^ since LO2 contributions can only affect the twist-stretch coupling 
of shape III and shape III already has twist-stretch couplings in both the 
isotropic and anisotropic models. 

We see from ( fL9|) that the twist-stretch coupling in our rigid plate model 
vanishes when the helix axis offset r vanishes. However, calculations for a 
pretwisted rod of material with a three-dimensional isotropic elasticity yield 
a twist-stretch coupling for zero helix axis offset provided the rod cross- 
section is not circular |18[ . To achieve consistency with our results, it seems 
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Table 2: Determination of the twist-stretch coupling for each shape of B- 
form DNA is shown below. The symbol TS signifies that there is a nonzero 
twist-stretch coupling. 



likely that the magnitude of this coupling decreases as the elastic constants 
coupling in-plane strains (i.e., u xx , u xy , u yy ) increase. We can thus view our 
calculation as the limit when these elastic constants diverge. 

We also see that in the isotropic limit (K 2 3 = 0) the twist-stretch coupling 
for shape III scales quadratically with the small helix axis offset as found 
previously in Refs. [§, |TJ]. However, when the K23 anisotropy is nonzero, 
the twist-stretch coupling scales linearly with the small helix axis offset. We 
note that the twist-stretch coupling in the isotropic limit vanishes for shape 
II since it has uj 2 = 0. However, shape III yields a twist-stretch coupling in 
both the isotropic and anisotropic cases since it has uj 2 7^ 0. This may be 
explained by considering the manner in which the plates are stacked in the 
two different geometries. In shape II the plates are stacked perpendicular to 
the molecular axis. In this case, the plates can move along the molecular axis 
without rotating in response to the stretching force. In shape III the plates 
are not perpendicular to the molecular axis, but are instead perpendicular 
to the local tangent dr/dn. In order to align the plates perpendicular to the 
molecular axis, they must twist. This implies that an isotropic rod bent into 
the shape of a helix will have a twist-stretch coupling, but carving a helical 
shape out of isotropic material will not produce a twist-stretch coupling. 



Shape 



Helix Axis Isotropic Model Anisotropic Model 
Offset (K 23 = 0) (K 23 0) 



I 

II 
III 
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7 Comparison to the Ribbon Model of DNA 



We compare our plate model to a ribbon model of DNA recently discussed 



in Ref. There we modeled DNA as a thin helical ribbon that is only 

allowed to stretch along the ribbon axis. Changes in extension along the 
molecular axis arose from changes in the angle the ribbon tangent makes 
with the molecular axis and from stretching the ribbon axis. We showed 
that the twist-stretch coupling scaled quadratically with the helix axis offset 
if the twist-bend, twist-stretch, and bend-stretch couplings for twists and 
bends about and stretches along the ribbon axis were set to zero. 

To make a direct comparison between the two models we consider the 
unstressed configuration 5 = (0, 0, 5 3 ) and to = (0,^2,^3) and do not allow 
stretch along the 1 and 2 axes of each plate by setting Ai = A2 = 0. This 
corresponds to a helical stack of plates with the planes of the plates per- 
pendicular to the tangent to the curve r connecting the geometric centers of 
the plates and stretching only allowed along the tangent dr/dn. When we 
set Ai = A 2 = in (]TT}) , we find a free energy identical to the one studied 
previously [T^]. We see below that our expression for the free energy ad- 



mits all of the couplings found in the ribbon model, namely the bend-stretch 
(f2 2 A 3 ), twist-bend (Q 2 Q 3 ), and twist-stretch (f2 3 A 3 ) couplings. (We have 
included the uj 2 contributions from the Jacobian and kept each coupling to 
lowest order in oj 2 /u) 3 and the second moments.) 



ribbon 

k B T = 

(A/28 3 ) J dnji^M + K 33 ( v l)Q 2 2 + ^ u ( % 2 ) + K 22 ( v l)^l + K m A 2 3 + 

2K 33 ^( V 2 )Q 2 A 3 - 2K 23 (r,l)n 2 Sl 3 - 2K 23 ^( V 2 1 )Q 3 A 3 \. (20) 
d 3 o 3 J 

We find that the presence of the the last two terms leads to twist-stretch 
couplings that scale linearly with the helix axis offset D ~ {K 23 / K 33 )rm/ S 3 . 
If we set i^23 = 0, the lowest order contribution is quadratic in the helix axis 
offset and we obtain the result D ~ -^33(^111/^3) 2 found previously |5|, [12| . 
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8 Conclusion 



We have presented an elasticity theory for B-form DNA modeled as a stack of 
thin, rigid plates. The symmetry of the DNA molecule and the assumption 
of plate rigidity enabled us to reduce the number of contributing elastic 
constants to four. We then derived an effective twist-stretch energy for DNA 
in terms of the relative excess twist a about the molecular axis and the stretch 
e along the molecular axis. We found an effective twist-stretch coupling with 
terms that scale linearly and quadratically with the small helix axis offset. 
We have argued that deformations of the base pairs are more costly than 
deviations in the shifts and rotations of the base pairs from their unstressed 
values. We therefore conclude that in the force regime where deformations 
of the plates are negligible, the twist-stretch coupling is due to a nonzero 
helix axis offset. We have also found that shape II, in which the normals to 
the base pairs parallel to the molecular axis, does not have a twist-stretch 
coupling in the limit K23 = 0. Since recent experiments show that there is 
a large twist response to stretch, this implies that either the coupling K 2 j, is 
relevant for an elastic description of B-DNA or that the unstressed straight 
state of B-DNA is composed of base pairs that tilt relative to the molecular 
axis. 

An important next step is to estimate the elastic constant K23 to deter- 
mine the relative magnitudes of the isotropic and anisotropic terms in the 
twist-stretch coupling. One might also calculate the twist-stretch coupling 
for DNA plasmids where the unstressed state of the molecular backbone is 
circular rather than linear. This calculation could provide a theoretical es- 
timate of the twist-stretch coupling found in experiments on complexes of 
Rec-A protein with DNA plasmidQ. 
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Appendix: Translational Invariance of the Plate Origin 



We comment here on the invariance of our model under the choice of the 
plate origin. This is an important feature since a twisted stack of plates can 
be described by an infinite number of one-dimensional curves depending on 
the choice of the reference point on each plate. We proceed by considering a 
change in the origin, or, in particular a change in the coordinate system such 
that x = (r)i,r)2,s) — > x = (771 , 772 , s)+R where R is an arbitrary vector. While 
this translation does not affect the relative rotation of two consecutive plates 
uj, it does affect the relative translation. Consider the action of g = (uj, 5) on 
x: 

x—^x + 5 — ujxx. (21) 

We now consider the transformation in (|2T|) shifted by R and compare it to 
the transformation g = (uj, 5) written in terms of the new coordinate system 
x. 



x + R = x + 5 — ujxx 



X 



R) + 5 -uj x 



X 



R 



+ R 



(22) 



In order for g and g to act the same way, we must choose uj = uj and 
5 = 5 — uj x R. We note that not only does (uj, 5) transform as above, but 
in addition, due to the linearity of the transformation, (Q, A) transforms the 
same way. As a result, we see that the derivatives of interest in (|8|) transform 

as 



d z Ui 



[A, 



z ijk 



Aj + SijkQj [Vk — Rk) + tijk^jRt 



/S 3 , (23) 



where R% is the projection of R onto the rji — f] 2 plane. Since changing the 
base point on each slab only requires components of R in that plane, we 
see that the last two terms in (^) cancel and the relevant derivatives are 
invariant with respect to the choice of the origin on the plate. Thus we are 
free to choose the most convenient origin for any calculation of interest. 
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Figure 1: Schematic cross section of DNA taken perpendicular to the sugar- 
phosphate backbones, ei is perpendicular to the line connecting the two 
phosphate backbones and points from the geometric center of the plate to- 
ward the minor groove and e 2 points from the geometric center toward the 
backbone running in the positive e3 direction. Each slice has a coordinate 
system (x, y, z) with e x as the x-axis, e 2 as the y-axis, and e 3 as the z-axis. 
Also, each slice is labeled by an integer n. 
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